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Random Walks on deterministic Scale-Free networks: Exact results 
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We study the random walk problem on a class of deterministic Scale-Free networks displaying a 
degree sequence for hubs scaling as a power law with an exponent 7 = log 3/ log 2. We find exact 
results concerning different first-passage phenomena and, in particular, we calculate the probability 
of first return to the main hub. These results allow to derive the exact analytic expression for the 
mean time to first reach the main hub, whose leading behavior is given by r ~ y 1-1 /^ wnere y 
denotes the size of the structure, and the mean is over a set of starting points distributed uniformly 
over all the other sites of the graph. Interestingly, the process turns out to be particularly efficient. 
We also discuss the thermodynamic limit of the structure and some local topological properties. 

PACS numbers: 05.40.Fb,89.75.Hc 



I. INTRODUCTION 

In the last few years, one of the most studied topic in 
network theory has been the investigation of highly in- 
homogeneous structures. Networks with strong variabil- 
ity in local metric and topological properties have been 
shown to well represent many real structures, occurring 
in nature and in man-made systems. Condensed matter 
and soft materials often feature an inhomogeneous orga- 
nization in space [l| , and even engineered devices can be 
constructed to reproduce a highly varying arrangement, 
in order to obtain the designed physical properties Q. 
Such networks, being physically embedded in space, are 
finite dimensional and it is known that in this case in- 
homogeneous topology can strongly affects physical phe- 
nomena occurring on the network itself. 

Graphs are also used to describe the generic relation 
between a set of elements or agents, as it happens in 
complex networks theory in biology, social science, com- 
puter science and economy 0,01 . Then, these networks 
often feature infinite dimensionality, and inhomogeneity 
has been detected in several real systems 0, H, 0, U- 
One of the most studied structures in the literature are 
Scale-Free (SF) networks, which show a degree sequence 
scaling with a characteristic power law [3j, [9( . This im- 
plies that highly inhomogeneous regions can be present 
in the network. 

Random structures and stochastic approaches have ap- 
peared to be very useful in this framework, even if, in real 
cases, one always has to deal with a single realization of 
the disorder. Due to the strong variability of the topol- 
ogy in the samples, the quenched properties could not 
coincide with the behavior found in the typical cases, as 
described by the probability distribution. Therefore, the 
study of quenched samples and deterministic topologies 
is a very interesting task [H [U E 0, 0, d ■ 

Apart from the discussion on how and why structures 
with scaling degree sequences are so often encountered 
in nature [16[, they certainly represent a new and in- 
teresting class of inhomogeneous graphs. Original tech- 
niques have been developed in order to characterize these 



topologies and their effect on physical properties. How- 
ever, the general topological and metric features of SF 
graphs, on the global and on the local scale, are still not 
completely understood. 

Random walks are one of the simplest stochastic pro- 
cesses affected by the topology of a network and, at the 
same time, the basic model of diffusion phenomena and 
non-deterministic motion. They have been extensively 
studied for decades on regular lattices [T3|, fractal net- 
works and finite dimensional inhomogeneous structures 
[l8| , where they have been shown to be able to evidence 
a new and unexpected phenomenon arising in presence 
of strong inhomogeneity, namely the splitting between 
local and average properties [lj| . The richer topology of 
a generic, inhomogeneous and infinite dimensional graph 
can have a dramatic effect on the properties of random 
walks, especially when considering infinite graphs, which 
are introduced to describe macroscopic systems in the 
thermodynamic limit. Random walks represent not only 
a good model for diffusion phenomena on large complex 
networks, but also a direct way to characterize their large 
scale topological features, also in presence of strong in- 
homogeneity, and their influence on physical properties. 
Once these are known, they could also be used to fruit- 
fully design and engineer a network topology with given 
properties. 

In this paper we want to deepen the analysis by 
studying random walks on a specific Scale-Free topol- 
ogy, namely a deterministic Scale-Free network, built in 
a recursive way and featuring a scaling distribution for 
the hub degrees, with an exponent 7 = log 3/ log 2. For 
this deterministic structure, first introduced in [20(, some 
metrical and spectral properties have recently been in- 
vestigated in details [2lJ. On the other hand, how these 
properties are linked to diffusion processes on the net- 
work is still an open problem. 

Using the formalism of random walks generating func- 
tions 11711. we derive exact expressions for the first passage 
times [22j for a random walker starting from the maxi- 
mally connected hub and from the last generation "rims" 
of the deterministic SF network. In particular, we inves- 
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tigate their dependence on the size of the network and 
we derive an exact expression for the mean time to first 
reach the most connected hub, namely the mean time to 
absorption if we place a perfect trap on the most con- 
nected hub. This quantity displays an extremely slow 
behavior as the size of the network increases, given by 
r ~ ]/ 1-1 / 7 , where V denotes the size of the structure, 
and the mean is over a set of starting points distributed 
uniformly over all the other sites of the graph. Interest- 
ingly, the trapping process appears to be very efficient. 
We also obtain an implicit relation for the generating 
function of the first passage probabilities from the hub 
to rims, and viceversa, which provides some insight into 
the leading singularity for the generating function of the 
return probability to the hub. In particular, we prove the 
recurrence of the main hub by exploiting the connection 
between random walks and electric networks [23] . More- 
over, we compare the return probability to the main hub 
with the return probability to an end-node evidencing a 
very different asymptotic behavior. All the results are 
checked with extensive numerical calculations. 

The paper is organized as follows: In Sec. [IT] we give a 
brief mathematical description of the deterministic Scale- 
Free graph, introducing the language and the formalism 
we will use in the whole article. Then, in Sec. IHII we 
derive the recursive relation between the first passage 
probabilities from hubs to rims and vice versa; from these 
we calculate the exact expression for the mean time to 
first reach the most connected hub and its dependence 
on the volume of the network. In Sec. IIVI we calculate 
the return probability on the main hub and we discuss its 
recurrence in the thermodynamic limit. Finally, Sec. [V] 
is devoted to conclusion and discussion. 



II. DETERMINISTIC SCALE-FREE NETWORK 

A generic graph Q consists in a non-empty set V of 
nodes joined pairwise by a set of links C [24] . Here 
we consider a particular set of deterministic graphs 
{G g }g=o,i,2,..., first introduced in [2(J, which can be built 
recursively: At the g-th iteration one has the graph of 
generation g, denoted by Q g (see Fig. [1]). 

Starting from the so called root constituted by one sin- 
gle node labeled as i = 1, at the first iteration one intro- 
duces two more nodes i = 2, 3 and connects each of them 
to the root; the resulting chain of length three represents 
the graph of generation g = 1. We call B\ = {2,3} the 
set of sites added at the first generation and linked to the 
root. Then, at the second iteration one adds two chains 
of length three and connects each end node to the root, 
namely B-2 = {4,6,7,9}, so that the root will increase 
its coordination number from 2 to 6. Proceeding analo- 
gously, at the g-th iteration one introduces two replica of 
the existing graph of generation g — 1 and connects the 
root with each site making up B g . 

Hence, at the g-th generation the root turns out to be 
the main hub with coordination number 2(2 9 — 1), the set 
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FIG. 1: (Color on line) Iterative construction of the determin- 
istic scale- free network; the first four generations are depicted. 
Nodes belonging to the set B g are represented in brighter 
color. 



of all nodes has cardinality V g = \V g \ = 3 9 and \B g \ = 2 9 . 
As shown in [U, Q g exhibits (2/3)3 s ~ l "hubs" with de- 
gree 2 l+1 — 2, being i £ [l,g — 1]. As a result, the tail of 
the degree distribution is a power law P{k) ~ fc~ 7 with 
exponent 7 = log 3/ log2 ps 1.59. However, this does not 
hold for the so called "rims" contained in the sets B g . 
In this case one finds P(k) ~ (2/3) k = e _7fc , where 7 = 
log(3/2) « 0.405, which shows that the scaling nature of 
the rims is not scale- free but exponential [2l[ . The topo- 
logical and spectral properties of this graph have been 
deeply analyzed in [21] where, in particular, the average 
degree is shown to be (k) g = J2kP(k)/V g = 4 [!-( 2 /3) 9 ], 
namely it approaches 4 as g — - ► 00. Indeed, on the one 
hand the graph becomes very complex and there exists a 
set of few [o(l)] nodes whose coordination number grows 
indefinitely, on the other the number of nodes with degree 
< 2 become large indefinitely; as a result the average de- 
gree remains finite, conversely the second moment (k 2 ) is 
divergent. It should also be underlined that by increasing 
g the number of cycles grows fast; as a result, while the 
number of end-nodes grows linearly with the volume V 
(~ 3 s ), the average distance from the main hub increases 
only logarithmically with V (<~ g). 



III. GENERATING FUNCTIONS AND FIRST 
PASSAGE PROPERTIES 

The simple random walk (RW) on a graph Q is defined 
by the jumping probability pij between nearest neighbor 



3 



sites i and j: 



A, 



Pij 



where Zj 



ZiSij and z 7 ; 



(Z _1 A) 



J^ igV Aij is the coordination 



number of site i. Therefore, the probability of reaching 
in t steps site j starting from i is 



P(i,j;t) = ( P % 



(1) 



In the following we denote with P g (i,j]t) the proba- 
bility that a RW on Q g , starting from a site i reaches the 
site with label j at time t and F g (i,j;t) the probability 
that the same walk reaches the site j for the first time at 
time t. Moreover, we consider the probability B g (t) that, 
at the generation g, a walk starting from any site in B g 
first reaches the central hub and the probability H g (t) 
that a walk starting from the central hub first reaches 
any site in B g . 

Now, the following equations hold 



B g (t) = 5 -f + - g /2/2 H i^)B g {t ~ 1 - k) 



(2) 



and 



i=i k=i 



1=1 k=l 

Some comments are in order here. In a graph of genera- 
tion g, each rim is connected to g nodes of which one is 
the main hub; this accounts for the first term in Eq. [2] 
The remaining g — 1 links connect each rim to a minor 
hub from which one can reach the main hub only pass- 
ing through a rim; this explains the sum of convolutions 
in Eq. [2] Analogously, from the main hub 2 9 links out 
of 2(2 9 — 1) point directly to a rim; the remaining links 
connect the main hub to nodes corresponding to rims of 
graphs of generations I < g and from such nodes one can 
reach any site in B g only through the main hub itself. 

The generating functions corresponding to Eq. [2] and 
Eq. [3] read as 

oo g-1 

B g (z) = = - + ~ a £^(2)4(z) ( 4 ) 

t=o 9 ^ i=i 



H g (t) 



2'J 



and 



9-1 



H g (z) ^H g (t)z< = ^z + I^E 2 " 1 ^ 



t=o 



i=i 



(5) 

Notice that, as can be easily inferred from Fig.[TJ Bi(t) = 
Hi(t) — St.i, from which Bi{z) = Hi(z) = z. We also 
recall that, by definition, B g {l) is just the probability to 
ever reach the hub from any site in B g and, analogously, 
H g (l) is the probability to ever reach any site in B g from 



i = l. Of course, for finite g one has B g {\) = H g (l) = 1. 
With some algebra we rewrite Eq. 0] and Eq. [S] respec- 
tively as 



H g {z) 



B g {z) 



23-1 

2s- 1 z 



9-1 



J2^ 9 B(z) =1. 



(6) 



(7) 



i=i 



Moreover, by properly handling Eq. [5] and Eq. we can 
obtain the following two finite-difference equations cou- 
pled together: 



1 



1 



B g {z) 



B g+ i{z) 
2 



B g {z) 



= --H g (z) 



(8) 
(9) 



which can be combined together to get the rather sym- 
metric expression B g (z) |[B g+ i(z)] _1 — [_Bi(z)] _1 | = 

2H g (z){[H g+ i(z)}-'-[Hi(z)}- 1 }. 

It is interesting to notice that, since < H g (z) < 1, for 
any g and for any z, from Eq. [6] we have z/g < B g (z) < 
z/[g(l — z) + z], from which one gets B OQ (z) = 0. Hence, 
in the thermodynamic limit the probability to eventually 
reach the hub from B g is zero. 

From Eq. [6] and Eq. [7] it is possible to calculate recur- 



sively B g {z 
ations we g' 

B 2 (z 
B 3 (z 
B 4 (z 



:et 



and 



H 2 (z 
H 3 (z 
H 4 {z 



and H g (z); for instance, for the first gener- 

(10) 

(11) 
(12) 



2-z 2 ' 

z(3-z 2 ) 
9 - 8z 2 + z 4 



z(3-z 2 )(14- llz 2 + z 4 ) 
168 - 282z 2 + 145z 4 - 24z 6 + 



2 z 



3-z 2 ' 

4z(2- z 2 ) 
14- llz 2 + z 4 ' 

8z(2-z 2 )(9-8z 2 + z 4 ) 
270 - 435z 2 + 211z 4 - 31z 6 + z s 



(13) 
(14) 
(15) 



The generating functions B g (z) and H g (z) allow to 
calculate the average time t g taken by a random walk 
started on a site in B g to reach the hub and the average 
time taken by a random walk started on the main 
hub to reach any site in B g , respectively. In fact, for any 
arbitrary generation g we can write: 



(16) 
(17) 
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Hence, deriving Eq. [6] and Eq. [7] and recalling that for 
finite structures B g (l) = H g (l) = 1, we find the following 
coupled equations: 



s-i 



if = 2-2 1 - s + 2~ 9 ^2 i tf. 



(18) 
(19) 



Notice that the previous two equations could be obtained 
directly from Eq. [5] and Eq. [3] by multiplying both sides 
by t and summing over t — 0, oo. 
Now, from Eq.[H]and Eq.[H]we get 



(20) 
(21) 



l g+l 


-t B 


= 1- 


r 9 


2t H 




= 2- 


h C 9 



Such equations can be properly handled to get 



and 



ZL g+2 OL H 



3t g+1 — 3, 



whose solution is given by 



i-H 



4 /3 
3 \2 



and 



(22) 



(23) 



(24) 



(25) 



where we used t B = tf = 1. 

Interestingly, tf and tf both display the same expo- 
nential law: The average time to first reach the hub from 
any of the 2 9 sites making up B g and the average time 
to first reach any site in B g from the hub grows with the 
generation as ~ (3/2) 9 . Also, recalling that \B g \ = 2 9 
and that V g =3 g , we can write 2 9 = y lo s 2 / lo s 3 and get 

(3/2)5 = Vg/\Bg\ = V" 1 - 10 ^ 2 / 10 ^ 3 . 

The average times tf and tf found above are useful 
to calculate the mean time to absorption r g . In fact, let 
us define t g (i) the mean time necessary to first reach the 
hub from site i ^ 1, then r g = Y^ieV — -0 an< ^ 

we can write 



£ *jW = £ E *«w 



(26) 



»ev 9 \v 9 



namely 



r 9 = T 9-l- 



V -l Vn 



- e *»w> ( 2? ) 

iev 9 \v 9 _i 



where V' 9 = V g \ {1} and V g \ V s _i is the set of sites 
added at the g-th generation. We now notice that the 
last sum in Eq. [27|is just given by 



E 

iev s \v g -i 

9-2 

-E 

i=i 



tg(i) 



tf\Bg\ 



(t 



\B 



[r 9 (3'-l)+3^ 1 + 3'tf] 2°- l -\ 



(28) 



where the first two terms simply allow for \B g \ rims and 
\B g \/2 nodes bridging between two rims; then, from the 
remaining nodes making up V g \ V g —i one can reach the 
main hub by passing through any rim, possibly via a 
minor hub, and this yields to the last sum. 

By replacing in the previous equation the expressions 
for tf and tf found in Eqs. EH US we get 



, . 2 32 f9\ 9 8 on 

> t„{i) = h— - --3 9 - 

^ 9W 5 15 V2 ' * 

»ev g \v g _i 



-2' J 



IV? 3'-l 
n- 



i=i 



which, together with Eq. [27] yields 
(3 5+1 - 1)t 9+1 - 3(3 5 - 1)t 9 = - 



x 3 s 



The solution of this recursive equation is 



1 



39-1 



32 /9 X9 
~9~ U 



9 



(17 + 4 5 )3* 



whose leading behavior is given by 

' " a v a 1 



2< 
(29) 



(30) 



(31) 



with 1 — log 2/ log 3 w 0.37. Such an exponent is even 
lower than the exponent log 2/ log 3 = 1/(7' — 1) « 0.63 
found in [T(| for a deterministic scale-free network dis- 
playing the power degree 7' = 1 + log 3/ log 2 = 1 + 7. 
Indeed, the multi-fractal nature displayed by the graph 
under study gives rise to a remarkable efficiency for trap- 
ping on i = 1. 



A. Numerical calculations 

For a generic graph, given the corresponding adjacency 
matrix A and the coordination matrix Z, the numerical 
calculation of the mean time to absorption can be per- 
formed by exploiting a differential equation where the 
normalized discrete Laplacian A = AZ _1 — I appears 
ficl [TP . I25I [26[ . More precisely, for the topological struc- 
tures analyzed here, the Laplacian A g is a V g x V g matrix 
which depends on the generation g and we have 



i=2 



1. 



(32) 



5 



so that Eq. [35] simplifies into 




FIG. 2: (Color on line) Mean first-passage time t 3 for a simple 
random walker moving on the deterministic sale-free network 
as a function of the volume V g ; data from the exact analytic 
expression in Eq.|30](») are compared to the asymptotic form 
of Eq. [31] (dash-dotted line) and to the numerical solution of 
Eq. the continuous line is a guide to the eye). 



Therefore, the mean time to absorption averaged over all 
possible starting sites i ^ 1 reads as: 



V, 



i=2 9 i=2 i'=2 



(33) 



In Fig. [2] we compare the analytical results of Eq.l30land 
of Eq. [3T] with the numerical results obtained via Eq. [33J 



IV. RETURN PROBABILITY 

The results found in the previous section allow to 
deepen the analysis of the random walk problem on the 
determinist structure considered. In particular, the prob- 
ability F g (l, l;t) that a random walk started on the hub 
returns to the hub itself for the first time after time t has 
the form 

i=l g H i=l v > 

(34) 

and the pertaining generating function is 

oo g „i-l 

p g (z) ^ Yl F B 0-> !; *)** = z E 2^i^ i(z) - (35) 

t=0 i=l 



By replacing the expression for B g (z) appearing in Eq. [9] 
we get the telescopic sum J2i=i[^ l+1 / Hi+i(z)—2 l / Hi(z)], 



F g (z) = 2 



= 2 



2(2f - i) 
1 



23+ 1 



2f-l 



H g+1 (z) 
{29-l)H g+1 { z y 



(36) 



which highlights that the probability to first return on the 
hub at generation g directly depends on the probability 
to first reach any site in $ g +i starting from i = 1. We 
also notice that for finite g, F g (l) — 1, that is the hub 
is a recurrent point, as expected for any point on finite 
graphs 

Interestingly, from Eq. 34 we can explicitly derive the 
average time vz to first return to the main hub as a func- 
tion of t^: 



l g+l 



- 1 



1 - 2-3 



2(3 9 - 23) 
23-1 



(37) 



where in the last equality we used the result of Eq. [24l 
Therefore, for large structures we have t° ~ (3/2) 9 , 
which is the same leading behavior found for tg and t g . 
This result is also consistent with Kac formula 123 ac- 



cording to which t\ 
E; e v„ Z su/ Z 



[P(l,l;oo)]- 



2\C a 



-'811 



J S11- 



From F g (z) we can now calculate P g {z) = 
ESoP ff (l, !;*)«* - [l-Fgiz)]- 1 [H: using Eq.[5]and 



Eq. [36] we get 



Pgiz) 



(23-l)H g+1 (z) 



23[z - H g+1 (z)} 
From Eqs. I13|[T51 we can derive 



1 2 
2~ ^2' 

3 4(2 -z 2 ) 
4(l-z 2 )(6-z 2 )' 

15 8(2-z 2 )(9-8z 2 + z 4 ) 

16 (1 - z 2 )(126 - 109z 2 + 22z 4 - z 6 ) 



(38) 

(39) 
(40) 
.(41) 



We checked up to generation g 
written as 

P g (z) = (23 1) 



that P g (z) can be 



/(*) 



(l-z 2 ) ff (z)' 



(42) 



where /(z) and g(z) are even polynomial of degree 2 9 — 2 
(in fact all possible cycles have even length), both de- 
void of any factor (1 — z 2 ) and satisfying f(l)/g(l) = 
1/(39 _ 29). In Fig.[3]we show data for P 8 (l, 1; t) and for 
Ps(i,i]t), being i G V g \ V 9 _i and Zj = 1, both obtained 
numerically from Eq. [T] P g (i t r t t) = [(A g + In- 
terestingly, the decay displayed by the two probabilities 
considered is significantly different, which provides a fur- 
ther evidence of the strong inhomogeneity of the graph. 
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FIG. 3: (Color on line) Normalized return probability 
P g (i,i;t) — Pao(i,i;t) for a random walk on Q% and started 
at the main hub (i — 1) and at an end-node belonging to 
i £ V g \ Vg-i, as a function of time t; only even instants of 
time have been depicted. 



Once the thermodynamic limit is taken, a divergence 
in P g (z) for z — > 1 would imply that the main hub is re- 
current [19j . Recurrence can alternatively be proven by 
exploiting the connection between random walks an elec- 
tric networks. The escape probability can be determined 
by calculating the effective resistivity R e ff at fixed dis- 
tance from the hub, when all the links of the network 
are replaced by unit resistors and a unit voltage is ap- 
plied between the hub and the points at fixed distance 
from it. Then, taking the thermodynamic limit, the re- 
lation between the escape probability from the hub and 
the effective resistivity of the network reads [23| : 



lim — — , 

9 ^°° z en R eff{g) 



(43) 



where, we recall, Z gn = 2(2 9 — 1) is the coordina- 
tion number of the hub at generation g. Therefore, if 
Z gll R e ff(g) diverges in the thermodynamic limit, the 
hub is recurrent, while if it converges to a finite value, 
there is a finite probability of escape from the hub. On 
highly inhomogeneous graphs such as our scale-free net- 
works, however, a vanishing P esc may not necessarily im- 
ply that the mean number of visits to the hub is infinite 

In general, when the unit voltage is applied to the re- 
sulting circuit, by symmetry one can detect nodes which 
are at the same voltage (e.g. the rims) and short them 
together without affecting the distribution of currents in 
the branches, nor the effective resistivity. In this way, 
by applying the standard rules for the sum of resistors 
in scries and in parallel, one can build an equivalent net- 
work for the circuit, in terms of the total resistivity. To 
better exploit the symmetry properties of the network, 



we consider only even values of g, so that the points at a 
maximum distance from the hub are topologically equiv- 
alent (an analogous relation can be easily obtained if only 
odd values of g are considered) . With some algebra one 
then obtains: 



1 3/2-1 1 1 
Reff iff) = zl = wTT 



23/2 ^ 2 k 
k=0 



1 - 



1 



29/ 2 



(44) 



holding for even g. Hence, recalling that Z gn = 2(2 9 — 1) 
we derive that the main hub is recurrent. 



V. CONCLUSION AND DISCUSSION 

In this work we studied the random walk problem on 
a class of deterministic networks exhibiting both a scale- 
free, P(k) ~ fc -7 , and an exponential, P(k) ~ (2/3) fc , 
degree distribution. The latter holds for a subset of nodes 
called rims which are at a distance 1 from the main node. 
Adopting a generating function formalism we calculate 
the exact average time r to first reach the hub, where 
the mean is taken over all possible walks connecting a 
site % to the hub and over all starting nodes i. The lead- 
ing behaviour for r turns out to be r ~ y 1 - 1 /i ^ where 
V is the total number of sites making up the network. 
Analogous power-law behaviours were found previously 
for exactly decimable fractals |25y2fJ as well as for other 
deterministic scale- free networks [1(| and Apollonian net- 
works 

.11] 

, however the structure considered here results 
remarkably effective, being 1 - 1/7 « 0.37. The reason 
lays in the large number of cycles and in the multi-fractal 
nature of the the graph under study which determines a 
very short average distance from the main hub. 

In the second part of the work we focused on the prob- 
ability to first return to the main hub and we obtained a 
recursive formula for its generating function which allows 
to get some insight into the local topological properties 
of the hub. In particular, the recurrence of the hub is 
proved by mapping the graph into an electrical network 
whose effective resistivity provides the escape probabil- 
ity from the hub. Hence, in the thermodynamic limit, 
a splitting between average and local properties occurs: 
transience in the average is expected due to the diverging 
dimensionality, while the hub is locally recurrent. 

Finally, we compared the return probability to the 
main hub with the return probability to one of the outer 
nodes evidencing different asymptotic behaviors. This 
further highlights the strong inhomogeneity of the graph 
considered. 
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